KOLAKOSKI-(2m,2n) ARE LIMIT-PERIODIC MODEL SETS 



BERND SING 



Abstract. We consider (generalized) Kolakoski sequences on an alphabet with two even 
numbers. They can be related to a primitive substitution rule of constant length I, Using 
this connection, we prove that they have pure point dynamical and pure point diffractive 
spectrum, where we make use of the strong interplay between these two concepts. Since 
these sequences can then be described as model sets with ^-adic internal space, we add an 
approach to "visualize" such internal spaces. 



1. Introduction 

A one-sided infinite sequence u over the alphabet A = {1, 2} is called a (classical) Kolakoski 
sequence (named after W. Kolakoski who introduced it in 1965, see [12]), if it equals the 
sequence defined by its run lengths, e.g.: 

u = ,22, JA^ 2 ^1^ 22 s _l_, ,22, JA^ 2 ^Ll_, ... 
^ 2 2 1 1 2 1 2 2 1 2 ... = u. 

Here, a run is a maximal subword consisting of identical letters. The sequence u' = loo is the 
only other sequence which has this property. 

One way to obtain u of (1) is by starting with 2 as a seed and iterating the two substitutions 

1 i-> 2 . 1 i — ► 1 

a : 2 ^ 22 and a i: 2 ^ ^ 

alternatingly, i.e., a substitutes letters on even positions and <t 1 letters on odd positions (we 
begin counting at 0): 

2 ,_» 22 i-» 2211 i — ^ 221121 ■-> 221121221 i-» . . . 

Clearly, the iterates converge to the Kolakoski sequence u (in the obvious product topology), 
and u> is the unique (one-sided) fixed point of this iteration. 

One can generalize this by choosing a different alphabet A = {p, q} (we are only looking 
at alphabets with card(*4) = 2). Such a (generalized) Kolakoski sequence, which is also equal 
to the sequence of its run lengths, can be obtained by iterating the two substitutions 

, v q i > j q ^ q q 

(2) cr : „ and a 1 : 

pi — > p 1 pi — > q y 

alternatingly. Here, the starting letter of the sequence is p. We will call such a sequence a 
Kolakoski- (p, q) sequence, or Kol(p, q) for short. The classical Kolakoski sequence u> of (1) is 
therefore denoted by Kol(2, 1) (and to' by Kol(l, 2)). 

While little is known about the classical Kolakoski sequence (see [6]), and the same holds 
for all Ko\(p, q) with p odd and q even or vice versa (see [22]), the situation is more favorable if 
p and q are either both even or both odd. If both are odd, one can, in some cases, rewrite the 
substitution as a substitution of Pisot type (see [22, 3]), which can be described as (limit-) 
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aperiodic model sets. By this method, Kol(3, 1) is studied in [3] and shown to be a deformed 
model set. The case where both symbols are even will be studied below. 

It is the aim of this article to determine structure and order of the sequences Kol(2m, 2n). 
This will require two steps: First we establish an equivalent substitution of constant length 
for Kol(2m, 2n) and analyze it with methods known from the theory of dynamical systems. 
Then we conclude diffractive properties from this. 

Remark: Every Kol(p, q) can uniquely be extended to a bi-infinite (or two-sided) sequence. 
The one-sided sequence (to the right) is Kol(p, q) as explained above. The added part to the 
left is a reversed copy of Kol(q,p), e.g., in the case of the classical Kolakoski sequence of (1), 
this reads as 

. . . 11221221211221|22112122122112 

where "|" denotes the seamline between the one-sided sequences. Note that, if q = 1 (or 
p = 1), the bi-infinite sequence is mirror symmetric around the first position to the left 
(right) of the seamline. The bi-infinite sequence equals the sequence of its run lengths, if 
counting is begun at the seamline. Alternatively, one can get such a bi-infinite sequence by 
starting with q\p and applying the two substitutions to get a 1 (q)\a (p) in the first step and 
so forth. This also implies that Kol(p, q) and Kol(g, p) will have the same spectral properties, 
and it suffices to study one of them. 



2. KoL(2m, 2n) as Substitution of Constant Length 

If both letters are even numbers, i.e., p = 2m and q = 2n (with m / n, where we can 
concentrate on m > n by the above discussion), one can build blocks of two letters and 
obtain an (ordinary) substitution. Setting A = pp and B = qq, these substitutions and their 
substitution matrix M (sometimes called incidence matrix of the substitution) are given by 

A i — ^ A m B m , ^ r (mm 
(3) a : B ^ A n B n and M = { n n 

where the entry is the number of occurrences of j in o(i) (i,j £ {A, B}; sometimes the 
transposed matrix is used). A bi-infinite fixed point can be obtained as follows: 

B\A^ A n B n \A m B m ^ ... 

This corresponds to the unique bi-infinite Kol(2m, 2n) according to our above convention. 

A substitution g is primitive if the corresponding substitution matrix M is primitive, i.e., 
M k has positive entries only for some k £ N. Equivalently, g is primitive if there exists a 
positive integer k £ N such that every i £ A occurs in g k (j) for all j £ A. The vector £ with 
components i i = \g(i)\, for i £ A, is called the length of the substitution g. If all l i are equal, 
g is a substitution of constant length. For the substitution a of (3), we have 

2m 
2n 

which is therefore not of constant length (recall that m ^ n). 

We will also need some notions from the theory of dynamical systems, see [18] and 
[8, Chapters 1, 5 and 7] for details. Let g be a primitive substitution over A and u £ A z a 
bi-infinite fixed point of g (i.e., u = g k (u) for some k £ N). Denote by u k the fcth letter of u 
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(k G Z) and by T the shift map (i.e., (T(u)) k = u k+1 ). Let A be equipped with the discrete 
and A z with the corresponding product topology. If we set 



X(g) = {T*{u) \k£Z}, 

then (X(g),T) is a dynamical system. Since we require g to be primitive, this dynamical 
system is minimal (i.e., {T k (u) \ k G Z} is dense in X(^) for all u G (X(^),T)), does not 
depend on the chosen fixed point u (if more than one exists, which is possible in the two-sided 
situation) and has a unique probability measure \i associated with it. In other words, it is 
strictly ergodic. On the Hilbert space C 2 (X(g), fi), we have the unitary operator 

U: C 2 (X(g)^) - C 2 (X(g),p,), 
f -> foT. 

If Uf = e tX f for some / / G £ 2 (X(p), ju), we call e* A an eigenvalue of (X(^),T) and / the 
corresponding eigenf unction. The spectrum (of the dynamical system) is said to be a pure 
point dynamical spectrum (or discrete spectrum), if the eigenfunctions span C 2 (X(g), fi). If 1 is 
the only eigenvalue and the only eigenfunctions are the constants, the spectrum is continuous. 
It is also possible that it has pure point and continuous components. In that case it is called 
partially continuous. Two dynamical systems (X, T) and (Y, S) are isomorphic (or measure- 
theoretically isomorphic), if there exists an invertible measurable map <p : X — ► Y, almost 
everywhere defined, such that <p preserves the measure and the dynamics (i.e., (p oT = Soip). 

The spectral theory of primitive substitutions of constant length is well understood. By 
the following criterion, we know that the substitutions a of (3) are related to substitutions of 
constant length. 

Lemma 1. [5, Section V, Theorem 1] Let g be a substitution of nonconstant length t. If I is 
a right eigenvector of the corresponding substitution matrix M, then (X(g),T) is isomorphic 
to a substitution dynamical system generated by a substitution of constant length. □ 

Since the substitutions a of (3) fulfill 1 the requirements of this lemma, the next 
task is now to construct the corresponding substitutions of constant length. This is 
achieved by numbering the ^4's and B 7 s in (3), i.e., we make the substitutions A m B m — > 
A 1 ... A m B 1 ...B m , respectively A n B n -» A m+1 . . . A m+n B m+l . . . B m+n . Then, the former 
substitutions (3) induce 

j-^-j Ai—-A m B 1 ...B m i-> (^i"^ m Bi--B m ) m (^ m+ i-^ m+ „B m+ i-B m+ „) m 

A m+ i---A m+n B m+1 ...B m+n i ► (A 1 ...A m B 1 ...B m ) n (A m+1 ...A m+n B m+1 ...B m+n ) n . 



^Note that from (2), one can also construct a primitive substitution by distinguishing odd and even positions, 
e.g., for Kol(4, 2) we would get (we use ~ as mark for even positions) 

4 -> 4444 
4 -> 2222 
2 -> 44 
2 -> 22. 

Instead of (3), we would get a substitution with substitution matrix 

( m m \ 
0mm 
n n 
\ n n ) 

which is also primitive (M 2 has positive entries only), but does not fulfill the requirements of Lemma 1. The 
eigenvalues of this M are {0, 0, 0, m + n}. 
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From this we get substitutions of constant length m + n (the eigenvalue of the substitution 
matrix M in (3)) by parting the right sides in blocks of m + n letters. For example, let m = 2 
and n = 1. Then 

A X A 2 B X B 2 h-> A^Bj B 2 A X A 2 B X B 2 A 3 B 3 A 3 B 3 
A 3 B 3 h-> A l A 2 B l B 2 A 3 B 3 

and one extracts the following substitution of constant length 3: 

-^4^ I ^ _/l_2 
A 2 i ^ B 2 A 1 A 2 

(z\ B l ^ B 1 B 2 A 3 

[ } B 2 ^ B 3 A 3 B 3 

■A^ i ^ ./1_2 
B 3 i-> B 2 A 3 B 3 . 

In the same way, we get substitutions of constant length m + n from (4). Note that these 
substitutions are all primitive since (compare to (4)) in every block of 2m 2 + 2n successive 
letters (note that we use m > n) every letter of A = {A ± , . . . , A rn+n , B 1 ,..., B m+n } occurs, 
so if (m + n) k ° > 2m 2 + 2n, then M k ° has positive entries only (this holds for fco > 3). Note 
also that we can reduce the alphabet by one letter by identifying A 1 = A m+1 (A 1 = A 3 in 
the example (5)), because both A's always yield the same substitution. 

Let us now determine the positions of A 1 in the sequence u generated by (4). They are given 
by a • 2m + (3 ■ 2n for some a, (3 € Z (e.g., 0, 2m, 4m, . . . , 2m 2 , 2m 2 + 2n, 2m 2 + An, . . . 2m 2 + 
2nm, . . .). Therefore we get gcd{i | u i = u = A t } = gcd(2m, 2n) = 2gcd(m, n). The height 
h(g) of a primitive substitution g of constant length i which generates a sequence u is defined 
as 

(6) h(g) = max{fc > 1 | gcd(k,£) = 1 and k divides gcd{i | u { = u }}. 
Then the following lemma holds. 

Lemma 2. Let (X(g),T) be a dynamical system, where g is a primitive substitution of con- 
stant length £ and height h(g) . Then the pure point part of this dynamical system is isomorphic 
to the dynamical system (Z^ x Z/h(g)Z,r), where r is the addition of (1,1) on the Abelian 
group Z £ x Z/h(g)Z, i.e., the direct product of the l-adic integers 7L t and the cyclic group 
Z/h(g)Z of order h(g). Therefore the pure point dynamical spectrum is given by 

{e--+ 2 ^ fc,neZ,meN}. 

Note that 

(7) Z,~Z Pi x...xZ Pr , 

where p 1 ,...,p r are the distinct primes dividing I, see [16, Section 3.10]. 

Proof. The lemma is just a reformulation of [5, Theorem 11.13], compare with [18, Section 
VI.] and [8, Section 7.3]. □ 

Proposition 1. Suppose (X(o),T) has pure point dynamical spectrum, where a is the sub- 
stitution of (3). Then 

'Z m+n x Z/2Z,r) ifm + n is odd 
'Z m+n , f) if m + n is even, 

where r is the addition of (1, 1) and f the addition of 1. 



(X(a),T) 
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Proof. For the substitution (4) of constant length £ = m + n, we have already seen that 
gcd{i | u i = u } = 2gcd(m,n). Therefore, using (6), the height of this substitution is 2 if 
I is odd and 1 if £ is even. The dynamical system of the substitution (4) is isomorphic to 
(X(a),T) by Lemma 1, therefore they have the same spectrum. The remaining statement 
follows from Lemma 2. □ 

We want to show that the spectrum of a is indeed pure point. For this we use slightly 
different substitutions of constant length that we deduce from a. We substitute A m B m — > 
a x . . .a m and A n B n —^b 1 ...b n (so in (4) we build essentially blocks of two, e.g., a 1 = A 1 A 2 ). 
We get 

(os a 1 ...a m ^ (a 1 ...a m ) m (b 1 ...b n ) m 

[) b 1 ...b n ~ (a 1 ...a m r(b 1 ...b n r, 

which again gives substitutions of constant length t = m + n. They are all primitive sub- 
stitutions by the same argument as before (in every block of m 2 + n successive letters every 
letter occurs). In the case n > 1, we can reduce the alphabet by one letter by identifying 
a ± = b x I— > a x . . . a m a 1 . . . a n . So we have two cases, n = 1 with substitutions 

<h ^ a i a 2 a 3 ... a m _ 1 a m a 1 
a 2 i-> a 2 a 3 a 4 ... a m a x a 2 



(9) : { 



a rn-l ^ a rn~l a m a l ■■• a m-3 a m-2 a m 

a rn ^ a m &!&!... 6j b x b l 
b l ^ 0l a 2 a 3 ... a m _! a m 6 X 



and n > 1 with substitutions (it is cumbersome to write down such a 9 in general form, 
but we will investigate its structure in the next section). Now the height of 9 and 9 is 
always 1, because if n > 1 we get gcd{i | n = u { = a x } = gcd(m,n), and if n = 1 we get 
gcd{i | u = itj = a x } = gcd(m, m + 1) = 1. 

Let g be a primitive substitution of constant length ^ and height /i(^) = 1. One says that 
q admits a coincidence, if there exist a A; € N and j < l k such that g k (i)j is the same for all 
i £ A (the jth letter of each g k (i) is the same, i.e., g k admits a column of identical values). 

Lemma 3. [5, Section III, Theorem 7] Let (X(g),T) be a substitution dynamical system of 
constant length and height h(g) = 1. Then (X(g),T) has pure point dynamical spectrum if 
and only if g admits a coincidence. □ 

If a substitution has height h > 1 , one gets a substitution of height 1 by combining letters 
into blocks of h letters. If this new substitution has pure point dynamical spectrum, so has the 
original substitution of height h, see [5]. Obviously, we get the following: if the substitutions 
9 and 9 (which arise from (8)) admit coincidences, then the dynamical systems defined by a 
of (3) have pure point dynamical spectrum. 



3. Coincidences and Coincidence Matrix 

Let us first check 9 of (9) for coincidences. For this we begin by exploring the structure: 
9(a 1 ) has two a 1 's at position and m, 9{a 2 ) has an a l at position m — 1, etc. We get an 
in 9{a k ) at position m + 1 — k for 1 < k < m — 1. Similar arguments show that there is an a m 
in 9{a k ) at position m — k for 1 < k < m and at m — 1 in ^(b^. Now, 9(a m ) has 6 1 's at all 
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positions 1, . . . ,m. Furthermore, has a b 1 at position m and shares the first m letters 

with ^(a 1 ). Schematically, we get the following structure of 9: 



(10) 



a m-l ^ • a m a l 

a m ^ a m b 1 b 1 ... b 1 b 1 b 1 
b 1 i— > * * * ... * a rn b 1 



Here we have omitted (.) all letters that are not necessary and by * we denote the part that 
6{a l ) and 6{b]) share. We now check for pairwise coincidences, i.e., for ii,i2 € A we check 
whether there is a k € N and a j < £ k = (m + l) k such that a k {i\)j = a k {i2)j- 

So we pick h,i2 € A = {a 1 , . . . ,a m ,b 1 }, i\ ^ ii- Suppose i\ ^ a m (otherwise we inter- 
change i\ and 22)- Then i2 either equals a m or at least 6^2) has an a m (every 9(i), i € A has 
one). In the first case take k = 1, otherwise k = 2. Observe that there are m successive 6 1 's 
in 9{a m ). So, if wc look at 9 k {i\) and 9 k (i2), we get the following: On the one hand, there are 
m successive b^s somewhere in 9 k (i 2 ), say at positions j, . . . ,j + m — 1. On the other hand, 
in 9 k (ii), there is at one of these positions j, . . . , j + m — 1 either a b ± , and we have a pairwise 
coincidence, or an a 1 . Say there is an a 1 at j with j < j < j + m — 1. Then in # fe+1 (ii) and 
# fc+1 (i2) we have pairwise coincidences at positions j ■£,... ,j • £ + m (the *'s of (10)). 

From this pairwise coincidences we get a coincidence inductively: We start with two letters 
ii,i2 and after k\ < 3 substitutions we have a pairwise coincidence, say at j\. Now a third 
letter i 3 may have something else at 9 kl (i^)j i , but whatever it is, in 9 kl+k2 (/c 2 < 3) all 
three coincide somewhere at a position j 2 with j\ ■ £ k ' 2 < 22 < {ji + 1) ■ £ ki - Since there are 
card(^l) = m + 1 letters, we get a coincidence after at most 3 • m substitutions (i.e., there is 
a j < £ 3m such that all 3m (i) ■ are the same for all i <G A). 

The structure of 9 is different. We have A = {a 1 , . . . , a m , b 2 , ■ ■ ■ b n } and therefore 
card(^l) = m + n — 1. Let us first show an example, with m = 5 and n = 3: 



(11) 



Since the positions of two consecutive a 1 's in the sequence differ by at most m, there is an 
a 1 in every 9{i) with i E A (note that 9 is a substitution of constant length £ = m + n). 
Again we look for pairwise coincidences, so choose ii,«2 £ -4- Then there is (at least) one a 1 
in 9(i\), say at position ji, and (at least) one in 9(i2), say at position j'2. Since there can be 
more than one a 1 in either, we choose ji,j2 such that \j\ — J2I is minimal. We further choose 
i\,%2 such that j\ < 32 (in the case ji = j'2, e.g., i\ = a 1 and i2 = a 5 in the above example, 
we are already done). If we look at 9(i\) and 9(12), there are two cases each (and therefore 
four cases, if we look at the combinations): Either 0(ii) jl+1 = a 2 , . . . ,9(ii)j 2 = a j 2 +i-j 1 or 





-> a x 


a 2 


a 3 


a 4 


a 5 


a i 


a 2 


a 3 


a 2 h 


-> a 4 


°5 


°1 


a 2 


a 3 


a 4 


a 5 


a i 


a 3 h 


-> a 2 


a 3 


a 4 


a 5 




a 2 


a 3 


a 4 


a 4 h 


-> a 5 




&2 


&3 


a i 




^3 


a i 


a 5 h 


-» 6 2 


^3 


°1 


^2 




a i 


^2 


&3 


b 2 h 


-> a 4 


a 5 


°1 


a 2 


a 3 


a 4 


a 5 


Oj 


63 H 


-> 6 2 


^3 




&2 


h 




^2 


&3 
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= b 2> ■ ■ • .0(*l)j a = 6 j 2 +l-ji and either = Om+l+jx-ia- ■ ■ > ^fe^-l = a m or 

0(^2)^ = b m+1+ j 1 _j 2 , . . . , Ofo) j 2 _i = b m . This is all that can occur by the chosen minimality 
of J2 - Ji > 0. 

Now we examine the case where 0(ii)j 1+1 = a 2 , ■ ■ ■ , 9(ii)j 2 = a j 2 +i-j 1 and ^(^ 2 )ji = 
a m+i+ji-j 2 ' • • • ' ^(* 2 )j2-i = a m' ^ e wan ^ to show that 9 2 {ii) and # 2 (z2) have a pairwise 
coincidence. Let us look at the a{s in 9{a i ) only (we use again the above example, but the 
reasons given apply for arbitrary m, n): 



a l 


1 — > 


1 

a l 


a 2 


a 3 


a 4 


a 5 


2 

a i 


a 2 


a 3 


a 2 


1 — > 


a 4 


a 5 


3 

a l 


a 2 


a 3 


a 4 


°5 


4 

a i 


a 3 


1 — > 


a 2 


a 3 


a 4 


a 5 


5 

a i 


a 2 


a 3 


a 4 


a A 


1 — > 


a 5 


6 

a l 


&2 


^3 


7 

a i 


^2 


&3 


8 

a x 


a 5 


1 — > 


b 2 


h 


9 

a l 


&2 


h 


10 


&2 


^3 



First we number the a 1 's with 1, . . . ,2m (left to right in 9{a i ) and top (i = 1) to bottom 
(i = m)) and we will speak of the kth a 1 (with 1 < k < 2m) according to that number. We 
observe the following: 

• Let k < m. If the kth. a 1 occurs at position j > n in 9(a i ), then the (A; + l)-st a x 
occurs at position j — n in 9(a i+1 ). If the fcth a 1 occurs at position j < n, then the 
(k + l)-st a l occurs at j + m in the same 9{a i ). 

• Let k > m + 1. If the A;th a x occurs at position j > n, then the (/c — l)-st a x occurs 
at j — n in the same 9{a i ). If the feth a 1 occurs at position j < n in ^(aj, then the 
(k — l)-st a x occurs at position j + m in 0(a i _ 1 ). 2 

• The second and the (2m)-th a l occur at the same position m in 0(a 1 ), respectively 
9(a m ). With the previous two observations we get: the kth and the (2m + 2 — /c)-th 
a l occur at the same position for 1 < k < m. 

• The first and the (m + l)-st a 1 occur in 9(a i ) where there is at least one more a 1 . 
This is obvious for the first a i , for the (m + l)-st observe that if it occurs at position 
j < m, then there is also one at j + n, and if it occurs at position j > m, then there 
is also one at j — m. 

These observations are based on the facts that the length of the substitution is m + n and that 
the position of the a 1 's in the sequence are separated by m or n only. Now the fact that the a i 
always occur in ascending order (i.e., we have a 1 a 2 a 3 . . . and not a 3 a 1 a 2 ... or something else) 
together with the first two observations essentially gives us an algorithm, which always yields 
a pairwise coincidence in 9 2 (i\) and 9 2 (i 2 ). Let us explain it in our example (12): Suppose 
we have i\ = a 2 and i 2 = a 3 . Then we have j\ = 2 and j 2 = 4. The first step is always to 
reduce j 2 by one, so we have j' 2 = 3. We have f 2 7^ ji, but there is a second a 1 in 9(a m ) (a m 
occurs at position j 2 in i 2 \), so we can increment j\ by 1 and get j[ = 3. We have j[ = j' 2 , and 
9(9{a 2 ).,) = 9(a 2 ) and 9(9{a 3 ).,) = 9(a 5 ) both have an a 1 at position 2 (the third respectively 

the ninth a^). Therefore we get a pairwise coincidence in 9 2 {a 2 ) and 9 2 (a 3 ). This algorithm 
relies on the "contrary line break property". 



Notice the contrary behaviour of the first two observations in going to a different or staying in the same 
0(a i ) and the position of the corresponding a x . We call this the "contrary line break property". 
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The other three cases are mutatis mutandis the same, see the positions of the a 1 's in (11). 
So, starting with any two 21,22 £ -4 we get a pairwise coincidence in 9 2 {i\) and 9 2 {i2). 
Inductively like before, we get a coincidence after at most 2 • (m + n — 2) substitutions. 
Therefore we have established the following. 

Theorem 1. (X(a),T) with a of (3) has pure point dynamical spectrum. Also, the dynamical 
system of the substitutions of constant length as defined implicitly in (4) and 9,9 of (8) and 
(9) have pure point dynamical spectrum. □ 

Proposition 1'. We have 

( 13 ) r x / a \ T ) ~ J x Z / 2Z > T ) ifm + nis odd, 

I (Z m+n ,f) ifm + n is even, 

where r is t/ie addition of (1,1) anc? f t/ie addition of 1. □ 

Remark: Let p be a primitive substitution of constant length ^ and height 1 over the al- 
phabet A = {1, • • • ,r}. Then we can define the coincidence matrix C, which is a quadratic 
ir-(r + l) x ir-(r + 1) matrix. The entries are defined as follows (where t < s and v < u): 

r _ j \{j I Q(s)j = u A g(t)j =u}\ ifu = v 

^(st)(uv) I |{j I g( 8 ). = u A ^(t) j = «}| + |{j I 5(5)^. = V A ^(i) j = if u + V 

Note that the substitution matrix M is a submatrix of C, since M su = C^ ss ^ uu y Also, C 
has row sums t. With this definition, Lemma 3 reads as follows. 

Proposition 2. [18, Proposition X.l] 3 For (X(q),T) are equivalent: 

(i) (X{q),T) has pure point dynamical spectrum. 

(ii) t is a simple eigenvalue of the corresponding coincidence matrix C. □ 

Obviously, t is an eigenvalue of C (C/£ is a stochastic matrix with row sum 1). 
Now the above proof of Theorem 1 translates to the following statements for C: 

• For 9, the third power of the coincidence matrix, C 3 , has a column (C^ st ^ a ^ a ^) with 
nonzero entries only. 

• For 9, the square of the coincidence matrix, C 2 , has a column (C 2 st ^ a a ^) with 
nonzero entries only. 

Lemma 4. [18, Lemma X.3] Let B be a quadratic matrix with nonnegative integral entries 
and row sums i. If Bij > 1 for all i and a fixed j, then £ is a simple eigenvalue of B. □ 



3 Note, however, that we use a definition of C different from [18]. The coincidence matrix there has dimension 
r 2 x r 2 and has the form (with the proper enumeration of the pairs) 




while the one defined above has the form 

M 
R P + Q 

Here M, P, Q are quadratic matrices and M is the substitution matrix. The Proposition is true for both 
matrices, the proof is analogous. 
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This establishes the desired result that £ is a simple eigenvalue for the coincidence matrix 
C of 6 (0), since £ 2 (I 3 ) is a simple eigenvalue of C 2 (C 3 ). 

We end this section with an example. We take 9 for m = 2, n = 1 and therefore the 
substitution 





a 1 h 


-> a 1 


a 2 


°i 


(14) 


a 2 h 


-> a 2 










-> a x 


a 2 


61 



We get the coincidence matrix 



( 2 


1 











\ 





1 


2 











1 


1 


1 




















1 


1 


1 


1 


1 








1 





V° 





1 


1 





1 J 



( 4 

2 
3 
1 
3 

V 1 



Here, already C 2 has columns with positive entries only. 
{0,0,1,1,2,3}. 



2 \ 
4 

3 
12 2 2 
2 10 
2 2 12/ 

The eigenvalues of C are 



4. Model Sets and Diffraction 

A model set A(Q) (or cut-and-project set) in physical space W 1 is defined within the following 
general cut-and-project scheme, see [17, 1], 



7T 



x H — > H 



\ U / 

1-1 dense 

r 

where the internal space H is a locally compact Abelian group, and r C M. d x H is a lattice, 
i.e., a co-compact discrete subgroup of R d x H. The projection ir int (r) is assumed to be dense 
in internal space, and the projection ir into physical space has to be one-to-one on r. The 
model set A(Cl) is 

A(fi) = {tt(x) \ xer, vr int (x) eO) c R d , 

where the window $7 C H is a relatively compact set with nonempty interior. 

Let u be a bi-infinite sequence over A = {1, . . . , r} and ^ : A — > C, i 1— > be a (bounded) 
function which assigns to every letter a complex number (the scattering strength). Then the 
autocorrelation coefficients rj(z) are given by 

1 N 

rt(z) = lim — — - — - > v(u„) ■ v(u„, z ) 

n=~N 

provided the limits exist. We write S z for the Dirac measure at z, i.e., 5 z (f) = f(z) for / 
continuous. Then the correlation measure 7 of u is given by 

7 = ^2v(z)S z , 

zez 
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and the diffraction spectrum 4 is given by the Fourier transform 7 of this measure. If 7 is a 
sum of Dirac measures only, i.e., 7 = YlkeS^k ' with a countable set S (for any choice 
of complex numbers (cjjg^), then u is pure point diffractive, i.e., the diffraction spectrum 
consists of Bragg peaks only. Also, the d fc 's are the square of the absolute value of the 
corresponding Fourier-Bohr coefficient at k and therefore non-negative (real) numbers. If 
there is no Dirac measure in 7, except one at position 0, <5o, which is determined by the 
density of the structure only, then the diffraction spectrum is continuous. 

For substitutive systems, the diffraction spectrum and the spectrum of the corresponding 
dynamical system are closely related, see [7, 19, 20, 13, 14]. 

Proposition 3. [13, Corollary 1.] Let g be a primitive substitution of constant length £ with 
height 1 over A = {1, . . . , r}, where u is a fixed bi-infinite word of g. Define U i = {j € 
Z I Uj = i} for all i G A. We have Z = U . . . U U r , where U denotes disjoint union. Then 
the following are equivalent: 

(i) g admits a coincidence. 

(ii) The U i 's are model sets for 

^int 

R x Z £ — > Z £ 

U / U 

dense 

r = {(z, z)\z G Z} — ► z 



TT 



u \ 
1-1 

z « — 



(iii) The sequence u and the sets 5 U i are pure point diffractive. □ 

Note that properties "one-to-one" and "dense" are obvious, the interesting part is that 
there exist relatively compact windows (with respect to the ^-adic topology) with nonempty(!) 
interior. We will discuss this point in the next section. 

With this proposition, we know that 9 and 9 (as defined implicitly in (8)) generate se- 
quences which are pure point diffractive. But we got every letter i £ A for the appropriate 
alphabet for 9, 9 by building four-letter blocks in the substitution rule (2), e.g., in (14) 
we have a 1 = 4444, a 2 = 2222 and b l = 4422. Such a deterministic substitution rule 
(i.e., Kol(2m, In) is local derivable from the sequence generated by 9, respectively 9) does 
not change the nature of the diffraction spectrum, only the Fourier-Bohr coefficients. The 
diffraction spectrum of Kol(2m, 2n) can be calculated from the one generated by 9, 9 as fol- 
lows: For Kol(2m,2n) we only have scattering strengths c' 2m and c' 2n . If we therefore choose 
Cj (i £ A with respect to 9, 9) according to its four-letter-composition in {2m, 2n}, then the 
diffraction spectrum 7 of 9, 9 is also a diffraction spectrum of Kol(2m, 2n), where Kol(2m, In) 
is realized as an atom chain with atoms not on Z but on |Z (we get the diffraction spectrum 
of Kol(2m, 2n) realized on Z by a simple rescaling with the factor 4). For the example (14), 



So we think of u as an atomic chain, where there is an atom of type u n at position n with scattering 
strength v(u n ). We represent this atom as v(u n ) ■ 5 n and therefore get a (countable) sum of weighted Dirac 
measures with autocorrelation 7. 

5 By this we mean the special choice of v, where we set c { = 1 and Cj = for all j / i. 



KOLAKOSKI-(2m,2ri) ARE LIMIT-PERIODIC MODEL SETS 11 

this means that we choose 

c fli = <4 ■ (1 + e"— + e~ 2_ r + e"^— ) = 

c a 2 = 4 • (1 + + e" ~ + e~ A ~) = 

c bi = c' 4 - (1 + e"-) + c' 2 • (e" 2 - + e" 3 -) = (1 - OK - 4). 

So in this case, the diffraction spectrum of Kol(4, 2) is given by the one of U b only. All 
Kol(2m, 2n) are pure point diffractive. 

Lemma 5. For a sequence u = (p(v), where v is a bi-infinite fixed point of a primitive 
substitution and 4> '■ A v — > A* u is a morphism (where u (v) is a sequence over A u (A v )), the 
following statements are equivalent: 

(i) The dynamical system of u has pure point dynamical spectrum. 

(ii) u has pure point diffraction spectrum. 

Proof. This is a (weak) conclusion of [14, Theorem 3.2]. □ 

We therefore obtain the following. 

Theorem 2. All Kol(2m, 2ri) have pure point diffraction and pure point dynamical spectrum. 

□ 

Remarks: The dynamical system of Kol(2m, 2n) is isomorphic to 

[1 m+n x Z/4Z,r) if m + n= 1,3 (4) 
(15) (Z m+n x Z/2Z,r) ifm + n = 2 (4) 

( z m+n> f ) ifm + n = (4), 

where r is the addition of (1,1) and f the addition of 1. This can be seen by the fact that 
we get Kol(2m, In) from the one generated by a of (3) by the substitution A = pp, B = qq, 
which corresponds just to doubling each letter in the latter one (see the substitution matrices 
in (3) and Footnote 1). 

If we consider how to get from (4), respectively (8), to Kol(2m,2n) and compare this to 
Proposition 3, we get (compare with (15)): Kol(2m, 2n), respectively U 2m , U 2n are model 
sets for 



7T 



J m+n " m+n 



\ u / 

1-1 dense 
r = {{z, z, z mod ord(F))jz G Z} 



where 
(16) 



' Z/4Z if m + n = 1,3 (4) 
Z/2Z if m + n = 2 (4) 
{0} if m + n = (4). 

(As before, Z m+n ~ Z p x . . . x Z p ^, where p 1 , . . . ,p r are the distinct primes dividing m + n.) 
The diffraction spectrum calculated from this cut-and-project scheme is consistent with the 
previous one, since the Fourier-Bohr coefficients which arise in each Z m+n separately are 

weighted by factors 1, e ^ , e~ 2 ~ or e~ 3 ~ which depend on the element of the cyclic group 
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Z/4Z (similar for the case Z/2Z), compare with [2] — but this is just how we calculated the 
Cj's from c' 2m and c' 2n . 

Let us show how one calculates the diffraction spectrum of U b of (14) explicitly. This sub- 
stitution can be written in recursive equations for U a ^, U a ^ and U b by observing at which 

position in which substitution a certain letter occurs (e.g., b 1 occurs in 9(a 2 ) at positions 1 
and 2 and in 0(6 1 ) at position 2) 6 : 

U ai = 3U ai U 3U ai +2 U 3U h 
U a2 = 3U ai + l U 3U a2 U 3U bi + l 
\ = 3U a2 + l U 3U a2 + 2 U 3^+2, 

where rC/j + s = {r-z + s|zG[/j}. Iterating these equations and using U a U U a ^ UU b = Z, 
one gets: 

U bi = (9Z + 5) U (27Z + 17) U (27Z + 22) U (81Z + 53) 
U (81Z + 58) U (81Z + 64) U (81Z + 65) U ... 

The Fourier transform of each lattice coset w rZ+s = J2z& K-z+s ls eas y t° calculate: 

Every Fourier-Bohr coefficients of U b are then given by the sum of the Fourier-Bohr coef- 
ficients of the corresponding u5^~^ s - The structure of U b (similar for all Uj^ that occur for 
substitutions of constant length) as a union of a countable but infinite set of (periodic) lattice 
cosets r • Z + s gives rise to the name limit-periodic, see [9] . 

The support of the Bragg peaks of Kol(2m, 2n) is given by 
k 

kez, seN 
k 



4 ■ (m + n) s 



2 s ■ Pi 1 ■ ■ ■ ■ ■ P7 



k e Z, si, . . . , s r € N , e e {0, . . . , log 2 (ord(F))} 



where p 1 , . . . ,p r are the distinct primes dividing m + n and F is the cyclic group of (16). 
However, there need not be a Bragg peak on every point of the support, e.g., Kol(8,4) is 
equivalent to a substitution 8 of constant length t = m + n = 6, but the positions of a 
letter a i; b i are separated by multiples of 2: The support in this case is better described by 
{|r | k S Z, s G N } than by | A; G Z, s,r£ N }. 



6 Note that these recursive equations form an iterated function system (IFS) in 3-adic space, because multipli- 
cation by a factor of 3 is a contraction in the 3-adic topology. The closure of the windows in the 3-adic internal 
space is therefore given by the unique compact solution of this IFS by (a generalized version of) Hutchinson's 
theorem [11, Section 3.1(3)]. This method is well known for unimodular substitutions of Pisot-type, see [15], 
[3] and the vast literature about Rauzy fractals. Similar results apply for all primitive substitutions of constant 
length £ in i?-adic space. 
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5. Euclidean Models of £-adic Internal Spaces 

So far, we have talked in an "abstract" way about the £-adic internal space. Usually the 
discussion ends at this point, but we want to "visualize" this £-adic space. We hope that by 
doing this, we also gain some intuition for such spaces and the meaning of p-adic internal 
spaces for model sets. 

Recall that a p-adic integer can be written as a formal series t = J2i>o *i ' P % w ^ n integral 
coefficients t i satisfying < t i < p — 1 (Hensel expansion). For the following, we identify a 
p-adic integer t with the sequence (ij)j> of its coefficients. The set of all p-adic integers (a 
ring) is written as Z , while the field of p-adic numbers is written as Q p and can be seen as the 
set of all Laurent series X^i>7V *i 'V % with N € Z. There is a p-adic valuation v p : Q p \ {0} — > Z 
defined by v p (t) = min{i E Z | t i ^ 0}, which gives rise to the p-adic metric with |i| = p~ v p^ 
and |0L = (and Q_ and Z are the completions of Q and Z with respect to the p-adic 
metric). So with respect to the p-adic metric, two numbers in Z are close if their difference 
is divisible by a high power of p. Note that this is a non- Archimedean absolute value (i.e., 
\x + y\ p < max{|x| p , \y\ p } for all x,y £ Q p ) and we therefore get some "strange" properties: 
all triangles are isosceles, every point inside a ball B r (x) = {y \ \y — x\ p < r} is the center of 
this ball, all balls are open and closed, etc., see [10]. 

For Euclidean models (see [21, Section 1.2]) of Z p we only need to know the formal series 
t = ^j>o ■ P % ■ We can even show models for Z^, where we do not make use of (7). For this, 
we use the addressing scheme known for fractals, for example in the Sierpinsky gasket, see [4, 
Chapter IV.]: 




Now the interesting thing here is that each point in the Sierpinsky gasket has a unique address 
(at least if we do not take the usual connected Sierpinsky gasket but the totally disconnected 
version; this can be obtained by using a contraction factor less than | in the IFS for the 
Sierpinsky gasket). So each point in the Sierpinsky gasket corresponds to a sequence (tj)j> 
with elements < t i < 2 — this is just the Hensel expansion of the 3-adic integers. Similarly, 
the Cantor set is such a geometric encoding of the 2-adic integers. "Reasonable" geometric 
representations of Z £ in M. d are those, where the sets Kr XQ x \ = {t € Z^ 1 1 = x , . . . , t r = x r } 
of points starting with the same address are represented by objects of the same size for a fixed 
r G N. Therefore we get that in d-dimension, Z^ with d + 1 < I < (kissing number in M. d ) + 1 
can reasonably be represented, if we do not make use of (7). Note that we can represent Z 3 
either in M 2 or R. 

This geometric representation surely fails for some p-adic (or £-adic) properties (all triangles 
are isosceles, every point inside a ball is its center, etc.), but some are also "preserved": points 
which are close in the p-adic topology are also close in this geometric representation and the 
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i linn hi hi hi nun ii 

Figure 1. 3-adic model for the internal space of (14) in R 2 (above) and 
R (below, stretched for better representation). The colors correspond to a 1 
(black), a 2 (dark gray) and b 1 (light gray). 

representation as totally disconnected fractal corresponds to the totally disconnected field 
Q p , TL V and its geometric models are both compact sets. And balls in the p-adic topology 
correspond to scaled down copies of the whole fractal. 

We like to conclude this section with our example from (14). The 3-adic geometric models 
are given in Figure 1. Observe that, in the two-dimensional representation, the parts (ac- 
cording to our above addressing scheme for the Sierpinsky gasket) K^ Q2 y, ^{12} anc ^ ^{21} 
are colored by only one color. This corresponds to the fact that at positions 9Z + 6 are a 1 's, 
on 9Z + 7 are a 2 's and on 9Z + 5 are b^s only in the bi-infinite sequence. So, large patches 
of the same color in the geometric representations correspond to lattice cosets l T 7h + s with 
small r. A similar addressing scheme can be used for the one-dimensional representation (and 
in fact for all £-adic representations). 
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